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Topics Today 


■ Vector Field Visualization (Kao) 

- CFD/numerical flow visualization 

- Particle tracking algorithms 

- Instantaneous vs time-dependent 
visualization 

- Unsteady flow volumes 

■ Surface Flow Texture Visualization (Shen) 

- Steady and time-dependent flow 
textures 

• LIC.UFLIC, LEA, IBFV, and IBFVS 

- 2D and 3D seed placement algorithms 

• Energy-based methods 

• Evenly-spaced methods 

• Topology-based methods 
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Vector Field Visualization - Applications 

Vector fields are found in many real-world applications: 

■ Study of flow around an aircraft 

■ Blood flow in our heart chambers 

■ Ocean circulation models 

■ Severe weather predictions 

The vector fields from these various applications can be visually depicted using 
a number of techniques such as particle traces and advecting textures 
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Flow visualization of a miniaturized heart pump designed to increase 
blood circulation in heart failure patients (NASA/DeBakey heart pump) 


Defense 


Astronomy (SOFIA telescope) 


Medical Applications 


Aerospace 


Summer Condition 


Winter Condition 




• Simulated airflow over across a mountainous portion of Northern California under two contrasting 
meteorological conditions 

• Visualization of cloud texture overlay with streamlines 

• Streamline paths are computed based on the input wind field 

• 4th-order Runge-Kutta integration scheme is used to advect the particles 

• Animations depicting flow over this geographical location provide immediate assistance in decision support 
and crisis management 
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Vector Field Visualization - Problem 

■ A vector field V (p) is given for discrete points p where p lie in either a 2D or 
3D grid 

■ 2D vector field visualization is straightforward 

■ 3D vector field visualization is challenging due to 3D perspective 

■ Time-dependent flow visualization has additional challenges 

■ A vector field V (p,t) is given for discrete points p and at many time steps 

■ Time-dependent flow visualization is also known as unsteady flow 
visualization 


Visualization Techniques 


■ Geometry-based methods 

rendering primitives built from 
particle trajectories 

- ID: streamlines, pathlines, 
streaklines 

- Line variations: tubes, ribbons 

- 2D: stream surfaces 

- 3D: flow volumes 
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Visualization Techniques - II 


■ Texture-based methods 

shading every pixels/voxels in 

the visualization 

-LIC 



- Texture Splats 

- Chameleon 



Visualization Techniques - III 



Topological-based methods 

extracting features based on 
topology 

- Critical points 

- Vortex cores 

- Skin-friction lines 
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Computational Fluid Dynamics 


A computational technology that 
enables one to study the dynamics of 
things that flow 

Runs on computer systems to model 
fluid flows using mathematically 
modeling, numerical methods, and 
software tools for pre-processing and 
post processing 


The wind tunnel has been traditionally 
used to simulate physical flows 


With the advent of more powerful 
computers, CFD is now the preferred 
means of flow simulation before final 
experimental testing (if any) 



A Space Shuttle orbiter model in NASA Ames 
Research Center’s 40x80 ft wind tunnel. 

NASA Ames Image Library (Photo: A. Melliar, 1975) 


Typical Steps in a CFD Analysis: Pre-Processing 

• Geometry processing (CAD surface geometry, surface domain 
decomposition) 

• Grid generation (surface grids, hole cutting, volume grids) 

• Flow solver input preparation (solver input files) 

• Boundary conditions (per grid basis, at each grid face, 
periodic, in-flow, out-flow, no-slip/viscous) 

• Flow conditions (e.g., Mach number, Reynolds number, angle of 
attack) 

• Numerical methods (parameters for the Navier-Stokes Eqs) 
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Typical Steps in a CFD Analysis: Flow Computation 

• Selection of physics model 

- inviscid/viscous/turbulence 

- body dynamics 

- no. of species 

• Selection of numerical methods 
parameters 

- accuracy 

- stability 

- convergence 

A corner view of the 5 1 ,200 core SGI Altix ICE system 
housed at the NASA Advanced Supercomputing facility. 

• Run the simulstion Photo Credit NASA Ames Research CerrterlMarco Libero 



Typical Steps in a CFD Analysis: Post-Processing 



As computation power continues to increase over the past decade, flow 
simulations require less amount of time to run. 

Nowadays, the pre-processing and post-processing stages consume most of 
the analysis time. 
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Common Methods in Experimental Flow 
Visualization - Adding Foreign Material 


• Streaklines: dye injected from a fixed position. By injecting the dye 
for a period of time, a line of dye in the fluid is visible 

• Timelines: a row of small particles (hydrogen bubbles) released at 
right angles to flow. The motion of the particles shows the fluid 
behavior 

• Pathlines: small particles (magnesium powder in liquid; oil drops 
in gas) are added to the fluid. Velocity is measured by 
photographing the motion of the particles with a known exposure 
time 



Experimental Flow Visualization 


NASA Photos (in Public Domain) 
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Experimental Flow Visualization 



Smoke and laser lighting sheet Oil flow visualization 

(NASA Langley, FS-2001 -04-64-LaRC) (NASA Dryden, IS-97/08-DFRC-02) 


NASA Photos (in Public Domain) 


Numerical Flow Visualization 



Though numerical flow visualization is not able to totally replicate the 
results from experimental flow visualization, it has been widely accepted 
as an effective mean to obtain accurate representation of the CFD flow 
solutions. 
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Numerical Flow Visualization - Basic Methods 
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Streamlines 


* Streamline: a line that is tangential to the instantaneous velocity 
direction (velocity is a vector, and it has a magnitude and a 
direction) 


• Release a particle into the flow and perform numerical integration to 
compute the path of the particle 
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Pathlines 

• A pathline shows the trajectory of a single particle released from a 
fixed location (seed point) 

• Experimental method: marking a fluid particle and taking a time 
exposure photo of its motion will generate a pathline 

• This is similar to what you see when you take a long-exposure 
photograph of car lights on a freeway at night 


Streaklines 

• Streakline: a line joining the positions, at an instant in time, of all 
particles that were previously released from a fixed location (seed 
point) 

• Continuously inject particles into the flow at each time step and 
track the paths of the particles 

• In a steady flow field, streamlines, pathlines, and streaklines are 
identical. However, they can be very different in unsteady flows 
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Timelines 

• Timeline : a line connecting a row of particles that released 
simultaneously 


• Timelines are generated by injecting rows of particles at some fixed 
time interval 




13 





2/16/10 


Stream/Path/Streak/Time Lines - Summary 

• Streamline : a field line tangent to the velocity field at an instant in 
time 

• Pathline\ shows the trajectory of a single weightless particle 
released from a seed point 

• Streakline-, a line joining the positions, at an instant in time, of all 
particles released from a seed point 

• Timeline-, a line connecting a row of particles that released 
simultaneously 

In a steady flow field, streamlines, pathlines, and streaklines are 
identical 


Particle Tracing 

The path of a massless particle at position p at time t is determined by: 

dp/c/t = v( p(t) ) for steady flow 
or 

dp/dt = v( p(t), t ) for unsteady flow 

The actual particle position at any given time can be determined by 
integrating the above equation (assuming that p(0) is given): 

t+At 

p{t + At) = p(t) + J* v(p(t),t)dt 
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Numerical Integration Methods 

• Euler’s Method: (Not recommended due to inaccuracy) 

Pk+i = Pk + v (Pk) Af p 0 = seed point 

• 2nd Order Runge-Kutta Method: (commonly used) 

P'=Pk+ v(Pk) Af 

Pk+i = Pk + MPk) + v(p) ]At/2 

• Higher Order Methods: (4th order RK and others) (recommended) 

a = 2Af vfpij, b = 2Af v(p k +a/2), c = 2At v(p k +b/2) 
d = 2Af v(p k +c/2), p k+1 = p k + (a+2b+2c+d)/6 

A good comparison of several integration methods can be found in [Teitzel et al. ’97] 


Notes on Particle Integration 

• The accuracy of particle tracing depends highly on the integration 
method used. 


• Flow solvers are often second-order accurate in time. Particle 
integration methods should be at least third-order or higher. 

• Though multi-stage integration methods such as 4th RK are 
commonly used, they require velocity data to be interpolated 
between the two consecutive time steps. 

• [Darmofal & Haimes ‘96] recommended using a fourth-order 
interpolant instead of linear interpolation: 

u"* U2 (x)= — u n +\x) + — u n (x)~— u nl (x ) + — u n -\x ) 

16 16 16 16 
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Basic Particle Tracing Algorithm 

1. Specify the seed point p(0), t=0 

2. Perform cell search to locate the cell that contains p(t) 

3. Interpolate the velocity field to determine the velocity at p(t) 

4. Advance the particle at p(t) using either RK4 or a higher-order 
method 

Note: cell search and velocity interpolation are needed for the 
intermediate points p* during the integration 

5. Repeat from Step 2 until the particle leaves the flow field 


Computational VS Physical Space 

• Particle Tracing in Computational Space 

- Map the given curvilinear grid and the associated velocity field onto a 
uniform grid by computing and then transforming Jacobian matrices 

- Since a typical grid is large, transformation is performed on-the-fly 

- Main advantage: cell search is eliminated 

- Disadvantages: Particle traces may not be accurate because the 
transforming Jacobian matrices are only approximations 

• Particle Tracing in Physical Space 

- More accurate because the interpolation and integration steps are 
performed in Cartesian space 

- Velocity field transformation is not necessary 

- Cell search is more time consuming, especial in multi-block grids 
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Cell Search Based on Tetrahedral Decomposition 

• Decompose the current hexadhedral cell into five tetrahedra 

• Odd or even configuration based on the sum of the node’s indices (i,j,k) 

• Coordinate and velocity interpolation are based on natural coordinates 

• Have been shown to gain a speed up factor of up to six times faster! 



Alternate between two configurations to ensure 

continuity between cells [Kenwright & Lane ‘96] 
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Particle Tracing Issues 


• Adaptive integration step size 

- Variable integration step size, df, is based on the local velocity vector 

• Particle tracing in multi-block grids 

- Overlapping multi-block grids are common in complex grid geometries 

- Particle paths are likely to cross multiple grids (grid jumping is necessary) 


• Particle tracking in moving grids 

- Particle position is based on the current grid block and time step 

- The Jacobian matrix needs to consider the grid velocity 


x n x t x r 

y% y rj y$ y r 

Zg Zyj Zg Zf 

ooo t M -t, 


(x x ,y x ,z t ) = p tM (§,r7,£)-p' f (§,rj,£) 


See [Lane *93] for more discussions 


Particle Traces - A Comparison 

• 2D oscillating airfoil 

• Unsteady flow 

• Moving grid 

• 200 Time Steps 

• Data set: Sungho Ko (‘95) 
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Instantaneous Flow Visualization 

• Streamlines, vector plots, contours, and cutting planes 

• Calculation is based on one time step of the flow 

• Effective for depicting the flow at an instant in time 

• Interactive visualization is possible 

• Time-variable is not used in the calculation 


Time-Dependent Flow Visualization 

• Calculation is based on many time steps 

• Time variable is used in the calculation 

• Effective for depicting time-varying phenomena such as vortex 
shedding, vortex breakdown, and shock waves 

• Require large disk storage 

• Interactive visualization is limited 
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Animation 


• Animation is an effective method to depict time-varying 
phenomena in unsteady flows 

• Sometimes, instantaneous visualization techniques are applied to 
unsteady flow data one time step at a time, then the results are 
animated 

• Instantaneous and time-dependent flow visualization can reveal 
very different information 

• Time-dependent flow visualization should be used to complement 
instantaneous visualization 
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Some Extensions of 3D Particle Tracking 

Stream surface [Hultquist ‘92] 

Stream ribbon and tube [Ueng, Sikorski and Ma ’96] 

Stream ball and streak ball [Brill et al. ‘94] 

Dash tube [Fuhrmann and Groller ’98] 

Time surface [Westermann, Johnson and Ertl ’01] 

Smoke surface [von Funck et al. ‘08] 

Integral surface [Garth et al. ‘08] 

Streak surface [Burger et al. ’09] and [Krishnan et al. ‘09] 

Flow volume [Max, Becker and Crawfis ‘93] 

Unsteady flow volume [Becker, Lane, and Max ’95] 


Flow Volume 

• Advect multiple streamlines from a polygon 

• Flow volume is the volume bounded by the 3D streamlines 

• Sorting is avoided by assuming monochrome flow volume 

• When an edge is too long, it is subdivided at midpoint 

• The flow volume is decomposed into tetrahedron cells 



[Max, Becker and Crawfis ‘93] 
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Unsteady Flow Volume 

• Advect multiple streaklines from a polygon 

• Flow volume is the volume bounded by the 3D streaklines 

• Flow volume evolves in time by moving “sideways” as the 
streaklines shift with the time-varying flow 

• All vertices in the flow volume need to be advected over time 

• Unsteady flow volume changes more rapidly than steady flow 
volume 


Adaptive Subdivision 



• Unlike steady flow volume, where only the last layer of particles (at 
the flow font) are advected in pseudo time 

• Unsteady flow volume requires the ENTIRE flow volume to be 
reconstructed at each time step 

• A new particle is inserted at the midpoint of a segment when a 
segment becomes too long 

• An existing particle is deleted when two adjacent segments become 
too short 
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Unsteady Flow Volume Over A Wing With 
Oscillating Flap 



[Becker, Lane and Max ‘95] 


Unsteady Flow Visualization/Large-Scale Data 
Visualization 

• Unsteady flow visualization often has the challenges associated 
with large-scale data visualization due to the time-varying nature of 
the data 

• Some previous workshops/symposiums on large-scale data 
visualization: 

o Symposium on Visualizing Time Varying Data (ICASE/NASA 
LaRC ‘95) 

o Visualizing Large-scale Data, BOF (Visualization '97) 

o Visualizing Large Datasets: Challenges and Opportunities 
(SIGGRAPH '99) 

o Time-Varying Data Visualization Workshop (Supercomputing ‘05) 

o Ultra-Scale Visualization Workshop (SC ‘06 - SC ‘09) 
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Large-scale Data Visualization Challenges 


• ‘Time-varying datasets present difficult problems for both analysis 
and visualization. For example, the data may be terabytes in size, 
distributed across mass storage systems at several sites, with 
time scales ranging from femtoseconds to eons.” 

'VDavid Banks et at., Symposium on Visualizing Time-Varying Data at 
ICASE/NASA LaRC '95 


• ‘The output from leading-edge scientific simulations is so 
voluminous and complex that advanced visualization techniques 
are necessary to interpret the calculated results. Even though 
visualization technology has progressed significantly in recent 
years, we are barely capable of exploiting terascale data to its full 
extent, and petascale datasets are on the horizon.” 

❖ Kwan-Liu Ma, Ultra-Scale Visualization at SC’07 


Progress in Large-Scale Data Visualization 

Though, there has been many good progresses made in: 

■ Flow visualization techniques 

■ Hardware-assisted/GPU-based visualization 

■ Out-of-core algorithms 

■ Parallel algorithms 

■ Multiresolution and data compression techniques 

■ High performance computing 

■ In situ visualization and data reduction 

The size of the flow simulation continues to increase . . . 
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Grid Systems of Past and Present 


Mid 1990's: 5 - 10 grids, 5 - 10 million volume grid points 



Courtesy of William Chan (NASA) 
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